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Abstract 

Recently, there appeared many papers devoted to the modeling of X-ray properties of 
Cygnus X-l and other black hole accretion disk candidates: e.g. J.Poutanen, J.Krolik, F.Ryde, 
MNRAS 292 (1997) L21; E.Agol, J.Krolik, ApJ 507 (1998) 304; A.Beloborodov, in "High 
Energy Processes in Accreting Black Holes", ASP Conf. Series, 161, (1999), p. 295. The goal 
of this electronic publication is to draw attention to our old papers where many ideas of recent 
discussions were anticipated (hot coronae, Comptonization, photon damping of waves, particle 
acceleration and matter ejection from accretion disks with large scale poloidal magnetic fields, 
etc.) 
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Heat transfer in the region of maximum energy release of an accretion disk will take place mainly by 
convection, serving to enhance the turbulence and to generate a powerful acoustic flux. The hard X-rays 
emitted by Cyg X-l (E <^ 200 keV) might result from Comptonization of soft photons in a corona formed 
around the disk through this heating. 

PACS numbers: 98.60.Qs, 97.70.Ss 

With high probability, the X-rays emitted by the source Cygnus X-l come from a black hole 
that is in a regime of disk accretion. The theory of disk accretion has been developed by several 
authors [l]-[4]. But difficulties arise if one seeks to explain the spectrum of Cyg X-l on the basis 
of this theory. The source has a luminosity L <j 10 38 erg/sec « 0.1L C , where L c = 1.3 • 10 38 M/M Q 
is the critical Eddington luminosity. For this luminosity the theoretical spectrum should fall off 
sharply in the energy range hv ~ 7 keV, whereas the observed spectrum extends up to energies 
as high [5,6], as ps 200 keV, or even beyond [7,8]. 

To overcome this difficulty, models have been proposed [9,10] wherein the properties of the 
radiating regions are essentially determined by the observations. An attempt has been made [11] 
to construct a self-consistent model for Cyg X-l that would account for the hard part of the 
radiation. This model, which neglects radiation pressure from the very outset, appears to have an 
internal contradiction in that it implies a disk as thick as its radius and drift velocity comparable 
to the orbital velocity. 

Hard radiation up to 200-keV energy can be present in the spectrum of an accretion disk only 
in the event that regions with hot (T e « 10 9 K) electrons participate in forming the spectrum; 
harder photons would require still more energetic particles. In an analysis of accretion disks we 
have shown [12] that hot regions or coronae with T e « 10 9 K will form around them. 

The standard theory presupposes that radiative heat conduction is responsible for energy 
transfer in the vertical direction. In the region of maximum energy release, P ~ P r ^> P g , 
k = k cs , and the vertical structure can easily be determined analytically. The density in this 
region is found [3] to be independent of z, while the temperature declines toward higher z. Such 
a situation will evidently be convectively unstable, since the entropy will decrease in the vertical 
direction. The onset of convective instability will serve to equalize the entropy. The mean density 
in the disk will become approximately an order of magnitude higher than the density of the disk 
in the standard theory. 

The convective heat flow Q conv , the convection velocity v, and the excess AVT of the tem- 
perature gradient above the adiabatic gradient (we shall take the mixing length / to be the half- 
thickness zq of the disk) are given [13] by the equations 



Cygnus X-l 



Qconv = C p pv(z /2)AVT, 



(1) 



Q com (l + AP r /P g )GM 
2pC p T R 3 



nl/3 



V = 



2/3 
Z J 



(2) 



1 



AVT='"« 



2/3 



pC p 



1/3 



i?z - 5/3 . (3) 



_(l + AP r /P - g)GM_ 

Hence we readily find that the excess AVT in the temperature gradient amounts to no more 
than 20% of VT, while the convection velocity is v = 3 • 10 8 cm/sec for r = 10r 9 , a value close to 
the velocity v s of sound in this region. The heat flow is carried mainly by convection (Q ~ Q C onv); 
because of the high density the radiative flux is equal to « 20% of the total flux. 

In the convective disk the flow Q ac of acoustic energy in the vertical direction is given [14] by 

Qa,c - pv 3 (v/v s ) 5 ~ 10 21 ~ 22 erg • sec -1 • cm 2 (4) 

and is a quantity of the same order as the total energy flux. An amount of order (P g /P r )Q&c 
is expended in heating the outer layers to an optical depth r < 1. One other mechanism serves 
to heat these layers [12]. A particle in a region near the surface of the disk that is transparent 
to radiation will be subject to the influence of radiation from the whole disk, and not only to 
the local radiation pressure gradient. The force of the radiation pressures will accelerate particles 
in the transparent region above the disk. Calculations of the equations of motion for particles 
subject to radiative, centrifugal, and gravitational forces show that in the region with r < 1 the 
particles (protons and electrons) will acquire vertical velocities corresponding to the temperatures 
(1 — 8) • 10 8 K of protons for L ~ 0.1L C . Turbulent relaxation of the particle motions in the corona 
will tend to equalize the mean energies of the ions and electrons and to form a quasi-Maxwellian 
particle velocity distribution. Thus the combined action of the two heating mechanisms we have 
described will produce around the accretion disk a hot corona with T e ~ 10 9 K for L ~ 0.1L C . 
The existence of an analogous corona has been postulated phenomenologically by Price and Liang 
[10]. 

We now estimate the density p co of the corona and the amount of material it contains. Since 
the gas pressure varies continuously with transition from the photosphere to the corona, we readily 
find that at the base of the corona 



Pco ^ PsT s /T co ~ 10- 2 p s « 10- 5 g/cm 3 , (5) 

where p s , T s denote the density and temperature in the photosphere of the disk for the parameters 
of the source Cyg X-l in the region of maximum energy release. The surface density of the corona 
is determined by the condition r cs ~ 1 and is ~ 2 g/cm 2 , which is more than an order of magnitude 
lower than the density of the opaque disk. 

The hot corona readily affords an explanation of the peculiarities of the Cyg X-l spectrum. 
The soft X-rays at hv <; 7 keV, which comprise 70% of the total flux, are formed in the 
photosphere of the opaque disk. Some of the radiation (~ 10%) will pass through the hot corona, 
and Comptonization will generate hard radiation up to hv ~ 3kT e ~ 200 keV, which will amount 
to 30% of the total flux. In Cyg X-l the luminosity varies around L = 0.1L C . Under these 
conditions the region of the disk with P r 3> P g , where convective heating is important, will be 
spatially small and will convert ~ 10% of the soft radiation flux, in accord with the requirements 
that the observations impose on the model. 

The changes in the spectrum as the luminosity of Cyg X-l varies [6,15,16] exhibit the following 
characteristic behavior. As the total energy flux rises, the radiation in the soft part of the spectrum 
(hv <^ 7 keV) increases, but in the hard range (hv ^ 10 keV) it remains almost constant, or perhaps 
may even decrease slightly. We shall assume that the variations in the luminosity are associated 
with fluctuations in the power of accretion. As the mass flow M rises in the region with P r 3> P g 
the fraction of the acoustic flow (4) expended in heating the corona will decrease: 

(P g /Pr)Q ac ~ M-\ (6) 
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For L ps 0.1 L c , when acoustic heating predominates, the rise in M may cause some decrease 
in the heating of the corona and in the amount of hard radiation. At the same time the flux in 
the soft range is determined by the radiation of the disk photosphere and will increase with M. It 
is worth noting that in strong bursts of luminosity, when L reaches about 0.3L C , the heating will 
begin to be governed by radiation-pressure forces, so the temperature of the corona and thereby 
also the power of the hard radiation should increase along with the rise in M and in the total 
energy flux. 

The observational evidence for the spectrum at hv > 150 keV is considerably less reliable than 
for the soft range. For instance, Baker et al. [8] have measured the radiation to hv ps 10 MeV, 
but find an observed signal of only ps 1% of the background (they also report a deficiency of 7-ray 
photons in the direction of Cyg X-l). In the range up to 600 keV, Haymes and Harnden [7] find 
a break in the spectrum at hv ~ 150 keV, and the spectral index a changes from a = 1.9 at 
hv < 124 keV to a = 3.1 at hv > 154 keV. 

Even if a corona with T ~ 10 9 K is present in the disk accretion model, radiation cannot be 
formed at hv = 200 — 2000 keV. Electrons with T e = 10 10 K or fast nonthermal electrons would 
be needed for that purpose. There are in fact two ways to form such fast electrons. Both would 
require the presence of a magnetic field in the disk. A magnetic field could exist in the disk either 
through twisting of the lines of force by differential rotation [3,17,18], or through infall onto the 
black hole of magnetized material having a small angular momentum [19]. In the latter case a 
poloidal magnetic field would be generated. 

In the binary system containing the source Cyg X-l, some of the material flowing from the 
giant star is dispersed in space near the system. The attraction of the black hole will not only 
produce an accretion disk. A small proportion of material having a low angular momentum will 
fall into the black hole and be decelerated in the disk. If this deceleration takes place at radii of 
(10 — 30)r 9 and if a thin collisionless shock wave is formed (as is very likely in the presence of 
an azimuthal magnetic field) wherein the kinetic energy is transformed into thermal energy and 
T e ps Tj, then hot electrons with T = 10 10 - 10 12 K will appear. The inverse Compton mechanism 
of interaction of the disk radiation with these electrons can lead to the generation of hard radiation 
with hv ~ 200 — 2000 keV or even higher energies. 

Another mechanism for producing fast particles is analogous to the pulsar process. If mag- 
netized matter with low angular momentum falls into the black hole (in addition to the disk 
accretion), a strong poloidal magnetic field will arise [19]. By analogy to pulsars [20], rotation will 
generate an electric field of strength E ~ —(v/c)B in which electrons are accelerated to energies 
e w R(v/c)Be w 3 • 10 4 [5/(10 7 Gauss)] Mev where v/c « 0.1 and R w 10 7 cm is the characteristic 
scale. In a field B ps 10 7 Gauss, such electrons will generate synchrotron radiation with energies 
up to ps 10 5 keV. Just as in pulsars, it would be possible here for e + e~ pairs to be formed and to 
participate in the synchrotron radiation. 

The authors express their appreciation to Ya.B. Zel'dovich, A. F. Illarionov, and I. S. Shklovskii 
for helpful discussions and valuable comments. 
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The X-ray brightness fluctuations induced in Cyg X-l and active galaxy nuclei by convection and 
turbulence in the sub-photospheric layers are discussed in terms of the disk-accretion and supermassive- 
star models. The variability time scales should be comparable with those observed, but in the case of 
supermassive stars and disks it is difficult to obtain the observed amplitude of the fluctuations. 

PACS numbers: 98.70.Qy, 98.50.Rn, 97.10.Cv 

1. One noteworthy property of the X-ray source Cygnus X-l is the variability of its flux on 
time scales ranging from milliseconds to a fraction of a second [1]. This variability is generally 
attributed to the presence of a turbulent accretion disk in the Cyg X-l system. In particular, the 
flux variations have been explained by the rotation of a hot spot [2,1] or alternatively, by the onset 
of instability in a zone where radiation pressure predominates [3,5]. 

In this letter we shall consider another model for the variability. We shall outline more fully 
what properties of the flux variations are to be expected on the accretion disk model as a direct 
result of the presence of turbulence, and how convective instability should develop in a region of 
high radiation pressure [6,7]. Acoustic waves generated in the convection zone will escape into 
optically thin layers, and will not only induce variable soft X-rays in the photosphere, but will also 
be responsible for variable heating of the corona. Comptonization of the photospheric radiation 
by the hot electrons under conditions where both temperature and density are variable will lead 
to flux variations in the hard range, hv ^ 5 keV. This X-ray fluctuation mechanism, involving the 
emergence of waves into transparent layers, evidently is of the same wave nature as the variability 
of the ultraviolet excess in stars experiencing intensive convection, such as T Tauri and UV Ceti. 

2. The waves escaping into the transparent layers and producing variable radiation in the 
photosphere and corona will occupy a rather narrow frequency band. Physically, the reason for this 
circumstance is that media with a high radiation pressure and a nonuniform distribution of plasma 
along the z coordinate (across the disk) will serve as efficient filters, isolating a characteristic 
frequency range from the broader spectrum generated by convection and turbulence. Waves of 
low frequency and a wavelength exceeding the scale height of the atmosphere will not escape 
outside but will induce oscillations of the coronal atmosphere as a whole. On the other hand, 
under conditions where radiation pressure predominates, high-frequency waves will experience 
very severe damping because of radiative friction, and their role in heating the corona will be 
insignificant. High frequencies may, however, appear in the observations either through nonlinear 
transformation of low-frequency waves or due to inadequate processing of the observational data 
[8]- 

We have examined elsewhere [9] the propagation of waves through a medium with strong 
radiation pressure, followed by their escape into the atmosphere. Waves emerging into the 
transparent layers will have a phase and group velocity equal to the velocity of sound in gas, 
v g = (jPg/p) 1 / 2 = (jTZT) 1 / 2 , where 7 =5/3 or 7 = 1 according as scattering or absorption 
predominates, 1Z is the gas constant, and the temperature T of the equilibrium atmosphere is 
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approximately equal to the temperature T c of the photosphere. The characteristic frequency uj c 
of waves emerging into the atmosphere is given by the expression 



For the accretion disk model, the gravitational acceleration g at radius r is equal [10] to 

GM h 

where the characteristic thickness is 



and Ro is the radius of the inner edge of the disk. The quantity H/H c = hqH/qc represents the 
ratio of the radiation pressure force to the gravity; kq is the opacity, including both absorption 
and scattering [9]. In a spherically symmetric star of luminosity L and radius R we will have 
g = GM/R 2 , H = L/4irR 2 , and H/H c = L/L c , where L c = 4ttcGM/k is the Eddington limiting 
luminosity. 

3. To calculate the characteristic frequencies of the fluctuations for Cyg X-l, we shall adopt 
the convective accretion-disk models given in a previous paper [7] and by Shakura et al. [11]. For 
a black hole of mass M = IOMq, we find that in both models as the luminosity rises from 0.1L C 
to 0.3L C there is little change in the characteristic frequency oj c corresponding to the equilibrium 
temperature, as given by Eq. (1); the values of uj c are confined [9] to the range 5-40 msec. Note 
that if a corona with T c ~ 10 2 T e is present, waves whose frequency is to = uj c or even somewhat 
lower will be able to escape. The frequency range mentioned is in good accord with the observed 
time-scales of variability. Our mechanism can yield fluctuations weakly correlated in time, and 
can simulate the white noise derived from observational analysis [8] of the brightness fluctuations 



Despite the good agreement between the theoretical time characteristics and the observations, 
analysis of small oscillations in our model still does not suffice to explain the amplitudes of the 
variability we observe, because pulsations should, in general, occur independently in regions whose 
size is comparable with the velocity of sound multiplied by the pulsation period, or about one-tenth 
the diameter of the zone of maximum energy production. Small pulsations should accordingly be 
smoothed out. However, the mechanism we are proposing could operate in a highly nonlinear 
regime. Strong nonlinearity of the waves would be expected in the light of theoretical estimates 
[7], and is essential to the very existence of the corona. The frequencies obtained from linear 
analysis correspond to the characteristic growth times of strong nonlinear flares. This problem 
awaits further theoretical treatment. Indeed, we would point out that not even the observational 
situation is fully clear [8]. 

An evaluation of the proportion r] of the acoustic flux that emerges into the corona and 
produces heating shows [9] that as the luminosity rises by a factor of 3 from 0.1 L c to 0.3 L c , the 
value of 7] will drop by a factor of 5-10. Thus, in accord with our previous suggestion [6,7], as the 
luminosity increases the power of the corona will remain approximately the same, or may even 
decline somewhat. This situation will prevent the hard X-ray flux {hv > 5 keV) of Cyg X-l from 
changing appreciably as the total luminosity varies, in agreement with the observations [12]. 

4. Multicolor photometry of the nuclei of certain galaxies has revealed the presence of com- 
paratively short-period brightness fluctuations in the nucleus of the Seyfert galaxy [13,14] NGC 
4151 (~ 130 days) and in the BL Lacertae object [15] OJ 287 (~ 180 days). There is presently no 
consensus as to the nature of active galaxy and quasar nuclei. Three models have been explored 




(7) 




(9) 



of Cyg X-l. 
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[16,17]: a) a dense star cluster; b) a supermassive star; c) disk accretion onto a massive black hole. 
The variability mechanism proposed here can operate in the two models b and c, for in both cases 
the subphotospheric layers will be strongly convective, resulting in the formation of a corona with 
variable heating. Using Eq. (1) to estimate the characteristic period of the fluctuations along 
with the standard model for a supermassive star [18] we find that for a mass M = 1O 8 M0 and a 
radius 100i? g (R g = 2GM/c 2 is the gravitational radius), the period corresponding to uo c would 
be ss 160 days, a value consistent with observation. On the accretion disk model, such a period 
would prevail in the zone of maximum energy production [9] for a black hole of M = 10 8 Mq, if the 
luminosity L = 0.1L C . In galaxy nuclei the fast component typically exhibits an approximately 
constant period in conjunction with sharp changes in phase [15], a behavior which, it would seem, 
accords with a convective wave origin for such fluctuations. 

It is worth emphasizing, however, that both in the supermassive star model and in the model 
of an accretion disk around a supermassive black hole, the ratio of the radiation pressure to the gas 
pressure is far higher than in the case of accretion onto a black hole of stellar mass. Acoustic waves 
will therefore be damped much more strongly. Our calculations indicate [9] that in this event the 
emergent acoustic flux will comprise no more than 1% of the flux generated at large optical depth 
- well below the variability amplitude observed. The supermassive star model would be supported 
if strictly periodic brightness oscillations of constant phase were detected, associated with the 
rotation or with pulsations of such a star as a whole, as might be the case [14] for NGC 4151. 
But analysis of observations of several variable nuclei does not reveal any strict periodicities [19]. 
Most likely the observed variability should be modeled by a random process [20] (see, however, 
Ozernoi et al. [21]). The model of a dense star cluster with frequent supernova outbursts offers 
the best fit, in our opinion, to the irregular variability of galaxy nuclei. 
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